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Matching QCD and HQET at three loops 

A.G. Grozin^ 

'^Institut fiir Theoretische Teilchenphysik, Karlsruher Institut fiir Tcchnologie 

QCD/HQET matching for the heavy-quark field [T] and heavy-hght quark currents f21 with three-loop accuracy 
is discussed. 
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1. Heavy-quark field 

QCD problems with a single heavy quark Q can 
be treated in a simpler effective theory — HQET, 
if there exists a 4- velocity v such that the heavy- 
quark momentum is p — mv + k (m is the on-shell 
mass) and the characteristic residual momentum 
is small: k <^ m. QCD operators can be written 
as series in \/m via HQET operators; the coeffi- 
cients in these series arc determined by matching 
on-shell matrix elements in both theories. 

At the tree level, the heavy-quark field Q is re- 
lated to the corresponding HQET field Qy (satis- 
fying i,Q, = g„) by [SI] 



Q{x) 



2m 



Qv{x), 



D^ = I?^ - w^v • D . 



(1) 



The matrix elements of the bare fields between 
the on-shell quark with momentum p = mv -\- k 
and the vacuum in both theories are given by the 
on-shell wave-function renormalization constants: 






(2) 



(HQET renormalization constants are denoted by 
Z). The Dirac spinors are related by the Foldy- 
Wouthuysen transformation 



u{mv + k) = 



2m 



O 



e 



Uy{k) . 



Therefore, the bare fields are related by 



Qo{x) 



imv-x 
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Qvq{x) 



where the bare matching coefficient is 



zo 



yos/Am + l) (ni + i) 



7os/'„("i) „("i)\ 



(3) 



(4) 



(we use the covariant gauge: the gauge-fixing 
term in the Lagrangian is — (9^^Q'^)/(2ao), and 
the free gluon propagator is {—i/p'^){giiu ^ (1 ~ 
0'a)PtxPiy/p'^)j the number of flavours in QCD is 
Uf = ni + 1). The 0(1/to) matching coefficient 
in ^ is equal to the leading one, zq; this re- 
flexes the reparametrization invariance [5]. The 
MS renormalized fields are related by the formula 
similar to (jSj, with the renormalized decoupling 
coefficient 



z(/i) 



ZQ(a^+^'(/i),a("'+i)(/i)) 



Zq. 



(5) 



If there are no massive flavours except Q, 
then Zq — I because all loop corrections 
are scale-free. The QCD on-shell renormal- 
ization constant Zq contains the single scale 
m in this case; it has been calculated [5] up 
to three loops. The three-loop MS anomalous 
dimensions of Qv |6I7I and Q |8I9) are also 
known. We have to express all three quantities 
Z°4igl,"'+'\ai-'+% ZQ(ai"'+i)(A.),a("'+i)(M)), 
ZQ(ai"'^(^),a("''(^)) via the same variables, say. 



„(«') 



(m), 



,("i) 



(fi), see [TU]- The explicit result 



for the renormalized matching coefficient z(^) can 
be found in 1 . Gauge dependence first appears 
at three loops, as in Zq [6]. The requirement 
of finiteness of the renormalized matching coeffi- 
cient (III at e — > has allowed the authors of [6^ 
to extract Zq from their result for Zq. 

In the large-/3o limit (see Chapter 8 in [TT] for 
a pedagogical introduction) : 
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'ulfi, 



5(w) + 0(^) , (6) 



where /3 = ,3oas/(47r), 7 = 7oQ!s/(47r) H (dif- 
ferences of n;-flavour and (ni -\- l)-flavour quan- 
tities can be neglected at the l//3o order). The 
difference of the QCD and HQET anomalous di- 
mensions 7 — ^Q—^Q (it is gauge invariant at this 
order) and the Borel image S{u) are |12|13|11] 

7(/3) = -2fF(-/3,0) = 
Po 



2Cf 

S{u) 



P 



(l + /3)(l + |/3) 



A5(2 + /3,2 + /3)r(3 + /3)r(l-/3) 
F{0,u) ~ FiO.O) 
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This Borel image has infrared renormalon poles at 
each positive half-integer u and at m = 2. There- 
fore, the integral in ^ is not well defined. Com- 
paring its residue at the leading pole u = 1/2 
with the residue of the static-quark self-energy at 
its ultraviolet pole u = 1/2 [14], we can express 
the renormalon ambiguity of z{ii) as 



^z(p) 



3AA 
2~ 



(8) 



(A is the ground-state meson residual energy). 
The matching coefficient is gauge invariant at the 
order l//3o. Expanding 7(/3) and S{u) and inte- 
grating, we obtain confirm the contributions with 
the highest power of ni in each term in our three- 
loop result, and predict such a contribution at 
at. 



Numerically, in the Landau gauge at n; = 4 



z{m) = 1 



A af\m) 

3 TT 



(16.6629-4.5421) 






(153.4076 + 42.6271 - 61.5397) 
(1953.4013-t----) ^ ^ + 



^^_ias M„i2.1208^"^^'^"^^ 



134.4950 



' (4)/ ^' 




(1953.4013 



(9) 



(/3o is for ni = A fiavours). Naive nonabelian- 
ization |12| works rather well at two and three 
loops. Numerical convergence of the series is very 
poor; this is related to the infrared renormalon at 
u= 1/2. 

Now let us consider the relation between the 
MS renormalized electron field in QED and the 
Bloch-Nordsieck electron field. The bare match- 
ing coefficient zq = ZT is gauge invariant to all 
orders, see W. In the Bloch-Nordsieck model, due 
to exponentiation, \ogZ^ = (3 — a^^'')a'^^'' / (Ane) 
(where the 0-flavour a^°^ is equal to the on- 
shell a « 1/137). In the fuU QED, logZ^ = 
a'^^^'a'^^'' /{4:T:e) + (gauge- invariant higher terms) 
(see [T] for the proof; this has been demonstrated 
up to four loops by the direct calculation |15)). 
The gauge dependence cancels in log(Z^/Z^) be- 
cause of the QED decoupling relation a'^-'a^^-' = 
0,(0)0,(0)^ Therefore, the renormalized matching 
coefficient z(/i) in QED is gauge invariant to all 
orders. The three-loop result is presented in [I]. 



2. Heavy light currents 

Now we shall consider |2] MS renormalized 
heavy-light QCD quark currents 



J(m) ^ ^, ^(M)io , 3q = Qo^Qa : 



(10) 



where F is a Dirac matrix. They can be expressed 
via operators in HQET 



,V) = Cr(/i)j(Ai) + ^Y1 B^{^AO^{^l) 



J 



where 



^l;^!' 



j(m) == Z, ^(/i)jo , Jo = qo^QvO : 



(11) 



(12) 



and Oi are dimension-4 HQET operators with ap- 
propriate quantum numbers. The leading-order 
matching coefficients Cr have been calculated up 
to two loops |12|16j . 

There are 8 Dirac structures giving non- 
vanishing quark currents in 4 dimensions: 

r = i, ^, il, ili. (13) 

7i 7j. > 7i llV > 7i 7!l7_l , 7i 7±7lT > 

where 7" = 7" — ^i'"- The last four of them 
can be obtained from the first four by multiplying 
by the 't Hooft-Vcltman 7^^^. We are concerned 
with flavour non-singlet currents only, therefore, 
we may also use the anticommuting 75 ^ (there 
is no anomaly). The currents renormalized at a 
scale /i with different prescriptions for 75 are re- 
lated by [m 

(<il^''Q)^ = Zp[i,){q^f'Q)^, (14) 

{H^'l'^Q)^ = Za{^^) (qirrQ)^ , 

(?-75^St"7^iQ), == Zt(m)(9-75"^7["7''1Q),, 

where the finite renormalization constants Zp^a,t 
can be reconstructed from the differences of the 
anomalous dimensions of the currents. Multiply- 
ing r by 7^^ does not change the anomalous di- 
mension. In the case of F = ji°'jP\^ multiplying 
it by 7^^^ just permutes its components, and also 



does not change the anomalous dimension, there- 
fore, 



Zt(m) - 1 ; 



(15) 



Zp.AifJ-) are known up to three loops !17]. 

The anomalous dimension of the HQET cur- 
rent (fT2|) does not depend on the Dirac structure 
F. Therefore, there are no factors similar to Zp^A 
in HQET. Multiplying F by 7^'^ does not change 
the matching coefficient. Therefore, the match- 
ing coefficients for the currents in the second row 
of p^ can be obtained from those for the first 
row. In the v rest frame 



Zp{^J.) = 
ZA{^J.) = 



Ci(/x) 



ZrifJ-) 



C AC (^) 

'5 

C^Hv(/^) C^0^1^2^3(^) 

C^HV^o(/i) Cyl^2^3{ll) 

C^AC^S jp) ^ (7^3 (^) 

C^HV^3(/Z) CjO-y^J^il-L) 

C^AC^a^i(p) ^ C^o^i {^l) 

C^HV^O^l (/i) Cj2^3{ll) 

C^AC^2^3(At) ^ C^2^3(^) 
C^HV^2..y3 (/i) C^O^l(^) 



= 1. 



(16) 



In particular, C^^^{^) = C [a fs\{ij). In the fol- 
lowing we shall consider only the matching coef- 
ficients for the first 4 Dirac structures in (ITS)) . 

In order to find the coefficients Cr(/i), we 
equate matrix elements of the left- and right-hand 
side of ()lip from the heavy quark with momen- 
tum p — mv + A: to the light quark with momen- 
tum fc„: 



<q{kg)\j{^i)\Q{mv + k)> = 

c^{^i)<q{k,M^l)\Q,{k)> + o 



/v, /Cn 



(17) 



The on-shell matrix elements are 

<q{kq)\j{fj-)\Q{p)> = Ug{kg)r{p,kq)u{p) 

xZri(^)4/2zV2, 

<qikq)\jOj')\Qv{k)> = Uq{kq)t{k, kq)Uy{k) 

xZr\^,)z'J'zy\ (18) 



where T{p, kg) and r(fc, kg) are the bare vertex 
functions, and Zq differs from Zq because there 
are no Q loops in HQET. The difference between 
u{mv -\- k) and Uy(k) is of order k/m, and can be 
neglected. It is most convenient to use k = kq — 
0, then the 0(l/m) term is absent. The QCD 
vertex has two Dirac structures: 

r{m,v,0) =r -{A + B]!!). 

This leads to 

u{0)T{mv, 0)u{mv) = T{mv, 0) u{0)Tu{'mv) , 
f{mv,0) =A + B. 

The HQET vertex has just one Dirac structure. 
Therefore, 



Cr(/i) 



rimv,o)z-\t,)zl/'zy' 



1/2^1/2 



f(0,0)Z7i(^)Z, 



(19) 



If all flavours except Q are massless, all loop 
corrections to r(0, 0), Zq, and Zq contain no 
scale and hence vanish: r(0,0) = 1, Zq — 1, 
Zq — I. The quantities T{mv,0), Zq, and Zq 
contain a single scale m; Zq has been calculated 
up to 3 loops in [6], Zq in pO], and r(mf,0) in 
the present work [2 . The MS renormalization 
constants Zj [7] and Zj [.18i (for all F) are also 
known to 3 loops. 

If there is another massive flavour (c in 6-quark 
HQET), then r(0, 0), Zq, and Zq contain a sin- 
gle scale rric. The first two quantities have been 
calculated up to 3 loops in [12]; the last one is 
known from [10]. The quantities r(?Tii;,0), Zq, 
and Zq now contain 2 scales, and are non-trivial 
functions of a; = mc/m. The renormalization con- 
stant Zq has been calculated in this case, up to 
3 loops, in [20] (the master integrals appearing in 
this case are discussed in Ref. [21]). The other 
two quantities are found in this work [2]. 

The bare on-shell QCD quantities T{mv,0), 



Zq, and Zq are expressed via 5q 



(»/) 



(and m^7^ 



if it is non-zero; we re-express it via the on- 
shell mass Trie). They don't contain /i. The MS 
QCD renormalization constant Zj is expressed 

via as" (li). The bare on-shell HQET quanti- 



ties r(0, 0), Zq, and Zq are expressed via g^ 



(«/-i) 



and 



'-cO 



(they are trivial at nic = 0); we re- 



express m^Q^ via the on-shell mass nic (which 
is the same in both theories). These bare quan- 
tities also don't contain /i. Finally, the MS 
HQET renormalization constant Zj is expressed 



via as ' (^). We re-express all the quantities 

I dH]) viaai"f'^\^j.), see [10]. 
From equation of motion we have 



idaj" = idaJQ = TOojo = m{fi)j{fi) 



(20) 



where m{fj,) is the MS mass of the heavy quark 
Q. Taking the on-shell matrix element between 
the heavy quark with p = mv and the light quark 
with kq = and re-expressing both QCD matrix 
elements via the matrix element of the HQET 
current with F = 1, we obtain [121 



toC^, (/i) — m(/i)Ci(/i) 



(21) 



The ratio m{ii)/m has been calculated at three 
loops numerically [22] and then analytically [23] 
(the analytical results |23|6j were later indepen- 
dently confirmed in |24j , and then in several other 
papers); for rric ^ 0, m,(fi)/m has been found 
in [201. 

The matching coefficients have been calculated 
up to 2 loops in T2' , and to 3 loops in the present 
work [2 . Analytical expressions are long; numer- 
ically, at rric = and fi = m we have 



^^ - 

(2) ^ 

(3) _ 



7.55-1-1.09 = 8.64, 
-5.47-^3.06 = -2.41, 
-9.87 +1.53 = -8.34, 



C 



-14.13 -f 2.42 



41.70, 



64.74 - 



C[ 



Cf^ = -37.25 

C{f} = -88.92 - 46.34 



&^\ = -123.61 



75.34-38.16= 101.92, 
-10.72 + 29.74= -18.23, 
45.34= -89.92, 
63.57 + 63.22= -123.96 



(in the middle part of each formula, terms with 
descending powers of p)^^ are shown sepa- 
rately). Naive nonabelianization [H] works rea- 
sonably well. 
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At 771c ^ 0, results are expressed via the mas- 
ter integrals depending on a; = nic/ra [3T]. Their 
status is summarized in the Tables 1-4 in this 
paper. In the present work [2], we were able to 
obtain exact analytical expressions (via harmonic 
poly logarithms of x) for 0(1) terms in the mas- 
ter integrals 5.2, 5.2a, from the requirement of 
finiteness of the matching coefficients. Therefore, 
the Table 3 in "21^ should be now replaced with 
the following Table [T] (DE means the method of 
differential equations, and MB the Mellin-Barnes 
representation). Unfortunately, 0{£) terms in 4 
master integrals are still known only as truncated 
series in x (the entries x in the table). Therefore, 
the rric corrections to the 3-loop matching coeffi- 
cients are also known only as truncated series in 
X (or numerical approximations). 

We let's apply our results to the matrix ele- 
ments between a _B or B* meson with momentum 
p and the vacuum: 

<0| [qit^'Q)^ \B> = -imBlEifi) , 

<0| (g-a"^Q) IS*> - /|.(Ai)(p"e^ - /e") . 



(22) 



The corresponding HQET matrix elements in the 
V rest frame are 



<0| 



(975^^ 



Q.) Ii?(/c)> =-7F(/i) 



(23) 



where the single-meson states are normalized by 
the non-relativistic condition 

<Bik')\Bik)>^^=i27rfSik'-k). 



These two matrix elements are characterized by a 
single hadronic parameter F{^) due to the heavy- 
quark spin symmetry. From (|20p we have |12j 



Ib 



TUB 



m{fi) ' 



(24) 



where we may replace tub by the on-shell 6-quark 
mass ?7i, neglecting power corrections. 

Our main result is the ratio fB-'/fB- At nic = 



Ib 2 ^ ^ 



{Cftf + Cat A + Tprnri + Tprh) Cf 






{Cprpp + CfCaTfa + C\rAA + CpTpnirpi 
CpTprFh + CATpniTAi + CATpTAh 



+ TpTifrii + TpTiirih + Tprhh)Cf 



' (4)^ \' 
Oil (m) 



O at, 



A 



(25) 



where 



1 2i ^ 1^ 4 2 31 

rp = -IT log 2 Cs TT H , 

3 ^ 2^"^ 9 48 ' 

1 2, o 1^ 1 2 263 

TA = — TT log2 + -Ca + -n^ . 

6 *^ 4^-^ 6 144 



19 

ri = — . 
36 
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41 



9 36 



8 1 2 

rpF = ~ q'^* - g log^ 2 - -TT^ log^ 2 

43 ^4 43 2 289 
24^^ ~ 192 ' 

^7r2log2 2 



1080 
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24 
115 
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lOR *^ 48 ^'^ 12 
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817 4 2233 2 
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1 ^2. 899 
-C3-Y^C3 



-TT* 



12960 
16 

rAA = -^^i 



648 
2log4 2 + l7r2--'^' 



119 2, „ 5, 11 

7r2log2 + -C5 

54 ^ 6^ 144 

ooon 48^25 



1^ ^4 

TT 



/r2 log^ 2 

5. 11 ,, 343 
^^ 144 



'^.Fi 



3240 
16 
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-04 



27 
"^'- 9 
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-TT 



1728 5184 

-Alog42+± 
27 ^ 27 

25, 11 , 

TT 



^ -^^2 



TT-^ 



28 2 1 <. — . 
TT^ log 2 H C3 



.a4-±]^-- 
27 



11 4 , 179 2 815 
162 '^ ^ 864 



32 4 , 4„ 4 2, 2„ 

log*2 + -^2iog2 2 
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46 2i ^ r. 1 4 1439 2 

■ — TT^ log 2 + 5C3 TT TT^ 

27 ^ ^ 162 1080 36 
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TT^ log 2 C3 H 

27 ^ 18^ 3240 
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-TT 
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-a4 + ^log2 27 
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^ TT^ log^ 2 



17 422 

72 ^"sT 
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1 2 203 
"" ^ "27'' ~ 324 
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r-.. = -gC3 + ^^ +324 
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Numerically, 

B 



^ -4.40-1.97= -6.37, 



fUmt)-^^"^ 



Ib' 



-4.26 + 0.89^-3.37, 



B' 



Ib 

pT 
B- 



(3) 



= -51.67 - 42.21 + 16.33 = -77.55 . 



/|.(™.)^*'^ 



-38.53 



-34.69- 22.91 -M9. 07 



Naive nonabelianization [IT works reasonably 
well. 

Asymptotics of the perturbative coefficients for 
the matching coefficients at a large number of 
loops I ^ 1 have been investigated in Ref. [33] 
in a model-independent way. The results contain 
three unknown normalization constants A'0.1,2 ~ 
1. The asymptotics of the perturbative coeffi- 
cients for Jb'/Jb contain TVq and N2] in the case 
of m/ifh it contains only iVo: 



/b. 

fB 

2 

7 



(n+l) 



L=-5/3 



14 

27 



\ + o 



50 



-9/25 r 



1 + - 
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N2 
No 






(26) 



The coefficient of N2/N0 is about 0.08 at n = 2, 
and it seems reasonable to neglect this contri- 
bution. Neglecting also 1/n corrections, we ob- 
tain EHl 



/b- 
Ib 



(3) 



L=-5/3 



14 
'27 



•56.37: 



-29.23. 



(04 = Li4(l/2)). The result for fl,{m)/fB- is 
similar. 



Our exact result —37.787 agrees with this predic- 
tion reasonably well. However, 1/n corrections 
are large and tend to break this agreement. It is 
natural to expect that 1/n^ (and higher) correc- 
tions are also substantial at n = 2. 
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